New insight on pseudospin doublets in nuclei 



B. Desplanques^* S. Marcos^ 
1 LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, INPG, 

F-38026 Grenoble Cedex, France 
^ Departamento de Fisica Moderna, Universidad de Cantabria, 

E-39005 Santander, Spain 

April 2, 2009 



Abstract 

The relevance of the pseudospin symmetry in nuclei is considered. New insight is 
obtained from looking at the continuous transition from a model satisfying the spin 
symmetry to another one satisfying the pseudospin symmetry. This study suggests 
that there are models allowing no missing single-particle states in this transition, 
contrary to what is usually advocated. It rather points out to an association of 
pseudospin partners different from the one usually assumed, together with a strong 
violation of the corresponding symmetry. A comparison with results obtained from 
some relativistic approaches is made. 

PACS: 21.60.Cs Shell model, 24.80.+y Nuclear tests of fundamental interactions and sym- 
metries, 24.10.Jv relativistic models, 21.10.Pc Single-particle levels and strength functions 



Keyvi^ords: Pseudospin doublets; Atomic nuclei 



1 Introduction 



The concept of pseudospin symmetry (PSS) in the description of single-particle nuclear 
levels has been introduced 40 years ago [HE]. It was motivated by the observation that 
levels within the same major shell such as 2si/2 and 1^3/2 {N = 2), 2^3/2 and I/5/2 
(A^ = 3), etc, are rather close to each other and, in any case, closer than usual spin-orbit 
partners. Thus, the pseudospin symmetry could be more accurately satisfied than the spin 
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symmetry (SS). Approximate pseudospin symmetry has been also found and analyzed in 
deformed nuclei [3]- [5] before a general explanation were suggested and lengthily accepted. 

Some theoretical explanation supporting the existence of pseudospin doublets was given 
by Ginocchio [6], within the phenomenological relativistic framework originally proposed 
by Walecka for nuclear matter [7] and used later on for studying properties of finite nuclei 
[H [9]. Looking at the small components of the Dirac spinors [lOj, he noticed that the 
equation allowing one to determine them leads to degenerate pseudospin (PS) doublets 
in the limit where the sum of the scalar (S5) and vector (So) components of the nucleon 
self-energies vanishes. In this limit, the radial parts of the small components are identical. 
This result suggested him to propose that the approximate pseudospin symmetry is due to 
the fact that the sum of the scalar and vector self-energies, S5 and Sq, cancels for a large 
part [II]-[l5]. The attractive explanation of the pseudospin symmetry has motivated a 
considerable number of studies by other authors both to extend its result [IS]- [21] and 
to get a better understanding of its limitations [2S]-[Sn]- Among these last ones, it was 
noticed that the quasi cancellation of the self-energies in the sum + Eq was associated 
with an energy denominator which, actually, does cancel at the nuclear surface |25]-|31j. 
[M] . This indicates that some care is required in concluding from the cancellation of the 
E5 and So potentials alone [32] • Thus, the expected similarity of the radial part of the 
small components of the Dirac spinors was not strongly supported by the detailed analysis 
of various solutions [28t |3T] . It even turned out that the similarity of the radial part of the 
large components for spin-orbit partners was better satisfied while the conditions for the 
energy levels required by the symmetry in this case are apparently worse [211 |3T]. Last, 
but not least, the pseudospin symmetry supposes the existence of single-particle states 
that are not observed [2], i.e. the PS partners of the so-called intruder states (states with 
n = 1 and j = I + 1). For clarifying some of the above points, a model where the above 
mentioned energy denominator does not cancel would be useful. In this case, however, 
any realistic nuclear model would not admit bound states. 

In this work, we mainly concentrate on the possible missing states. For our purpose, 
we first work in the non-relativistic case and consider two interactions that evidence the 
spin and pseudospin symmetry. For simplicity, we assume that they are related to each 
other by a unitary transformation that commutes with the kinetic-energy operator. To 
study specifically the missing-state problem, we also introduce an interaction that is a 
superposition of the two above ones with relative weights, 1 — x and x. By varying 
the X parameter from to 1, one should see how energy levels in the spin-symmetry 
case are related (or unrelated) to those in the pseudospin-symmetry one and vice versa. 
Some extension to a relativistic approach is made with the idea of an application to 
clarifying Chen et al. results [22], which evidence a spectrum very similar to ours in 
the nonrelativistic approach. Tentatively, we also considered the case of finite potentials 
instead of harmonic-oscillator-type ones used in the above work. 

The plan of the paper is as followed. The second section is devoted to the general 
aspects of the non-relativistic interaction we are using for our study. Results for spectra 
obtained from two different interaction models are presented in the third section together 
with some discussion. Complementary information from examining wave functions is 
given in the fourth section. The fifth section is devoted to a comparison with results 
from relativistic approaches. Conclusions are given in the sixth section. Two appendices 
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contain some technical details. 



2 Interaction model 

The interactions with spin and pseudospin symmetry we start with may be written as: 

V, = Vo{r) , 

Vi = a-pVi{r) a-p, (1) 

where Vo(r) and Vi(r) are spin independent and only depend on the r variable (local 
potentials). The operator a-p appearing in Vi is an unitary one. It changes the states 
corresponding to the pseudospin-symmetry space (represented with the symbol") to those 
corresponding to the spin-symmetry space [10]. In this operation, the total angular mo- 
mentum, j, is preserved while the orbital angular momentum is modified by one unit, 
implying a change of parity. Moreover, the operator makes the interaction Vi non local. 
The corresponding Hamiltonian may therefore be written as: 

indicating that the spectra corresponding to the interactions a-pVi{r) a-p and Vi(r) will 
be the same, apart from the spin-orbital angular momentum assignments of course. Wave 
functions for the Hamiltonian Hi can be obtained from those for the Hamiltonian Hi = 

2 

^ + Vi(r) by applying to these ones the operator a-p. Apart from the angular-momentum 
factor, the solutions of Hi for spin-orbit partners are the same in both configuration and 
momentum spaces. As the operator a-p is local in momentum space, this equality will 
turn into an equality for momentum-space solutions of Hi corresponding to pseudospin 
partners. 

The total interaction that is useful for our purpose of looking at the transition of the 
spin-symmetry case to the pseudospin-symmetry one may be written as: 

V{r) = {l-x)Vo + xVi 

= {1 — x) Vo{r) + X a-pVi{r) a-p , (3) 

where x is supposed to vary from to 1. Actually, as the interaction Vi is spin dependent, 
we can recover a part of the usual spin-orbit splitting by taking slightly negative values 
of X. 

As mentioned in the introduction, we consider the particular case where Vi{r) = Vo(r). 
This does not diminish the relevance of results to be obtained here, making them more 
striking instead. Without performing any calculation, one can see that the states for the 
spin and pseudospin-symmetry cases, shown respectively in the l.h.s. and r.h.s. of Table 
[H are in a one-to-one correspondence. The n quantum number represents the order of 
the states and, apart from a possible conventional shift by one unit, it is not, necessarily, 
identical to the node number of the radial wave function. This relationship, which is 
often assumed implicitly in the literature, is specific of a local potential. The question of 
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spin symmetry 



pseudospin symmetry 



2pi/2, 2p3/2, 1/5/2, 1/7/2 ^ 2Si/2, 2J3/2, 14/2, 1^7/2 
2Si/2, 14/2, 14/2 2pi/2, lp3/2, 1/5/2 

lpi/2, lp3/2 ^ 15l/2, 14/2 

lSl/2 ^ lpi/2 



Table 1: Representation of spectra in the spin-symmetry and pseudospin-symmetry limit 
(l.h.s. and r.h.s. respectively): states on the left in each row correspond to standard 
major shells. For the interaction assumption Vi(r) = Vo{r) considered in this work, states 
on the r.h.s. in a given row have the same energy as those on the l.h.s. in the same row 
and the same total angular momentum (but a different parity). We anticipated that the 
quantum number n is preserved in making this correspondence. 

interest is whether there is a global continuity between states on the l.h.s. of the table 
and those on the r.h.s. when the interaction varies from the spin-symmetry limit {x = 0) 
to the pseudospin-symmetry one {x = 1). It is reminded that the current understanding 
of the pseudospin symmetry supposes that the pseudospin partners of the states on the 
very right in each row of the left part of Tabled] are absent, due to specific conditions of 
realistic nuclear potentials [14]. The argument relies on a mathematical analysis of the 
nodes of wave functions. The authors do not examine in what measure this absence could 
have a relation with the differences between realistic models and those satisfying exact 
PSS. As we shall see in this work, as one approaches the PSS limit, a new type of nodes 
could appear, which may affect their conclusion. 

It remains to specify the choice of the potential Vo(r)(= Vi(r)). We consider two choices 
inspired by the harmonic-oscillator and the Woods-Saxon potentials. The comparison of 
results for these two models evidences a little but important detail that could be relevant 
for the interpretation of pseudospin doublets. The absence of singular behavior for these 
potentials lets one to expect a continuous transition between states on the l.h.s. and r.h.s. 
of Table [H 



3 Results and discussion for the spectra 

We, successively, present in this section level spectra for the harmonic-oscillator (HO) and 
the Woods-Saxon (WS) models. Parameters are those appropriate to the description of 
a nucleus such as ^°Ca. This is followed by some comments about the spectra. 
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Harmonic-oscillator model 

Spin symmetry Pseudospin symmetry 




-0.2 0.2 0.4 0.6 0.8 1 



X 

Figure 1: Single-particle spectrum evidencing the spin- and pseudospin-symmetry cases 
and the continuous transition between them. The interaction is based on a harmonic- 
oscillator one: V = ^Muq (j^l — x) r'^ + x a-p (^-p)- The observed order of the 
spin-orbit partners is recovered for negative values of x. The figure includes on the l.h.s. 
two "standard pseudospin doublets", 1^3/2 and 2si/2 on the one hand, I/5/2 and 2^3/2 on 
the other hand. 

3.1 Harmonic-oscillator model 

For the potential V{r) inspired from the harmonic-oscillator form, the two terms Vo(r) 
and Vi(r) entering its definition assume the standard form: 

Vo{r) = V,{r) = ^Mu',r\ (4) 

As results for the energy levels only depend on the parameter ujq, it is sufficient to present 
them in units of this parameter. They are shown in Fig. [1] for x varying from -0.2 
to 1. Actually, calculations were done with a Gaussian potential deep enough to be 
approximated by the HO potential. The reason is that this is more convenient for our 
study, which is made in the momentum space due to the non-locality of the interaction 
model we are using when x 0. Moreover, the Fourier transform of the Gaussian potential 
can easily be performed. As can be observed from the figure, there is almost no departure 
with the results of the harmonic-oscillator model when a comparison is possible (at a; = 0). 
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While the spectra can be essentially given in terms of ujq, the value taken by this parameter 
is nevertheless relevant in obtaining wave functions that could be compared to those of 
the WS model. We use ujq = 7 MeV so that the spectra of the HO and WS models roughly 
overlap. 

Examination of the spectrum around x = —0.2 shows it a priori contains the two 
standard pseudospin doublets, 1^3/2, 2si/2 and I/5/2, 2^3/2, though the corresponding 
states are slightly apart in each case. The comparison with the spectrum at x = 1, which 
evidences an exact pseudospin symmetry, suggests a different pattern however. The states 
ld3/2 and 2si/2 (or I/5/2 and 2^3/2) tend to separate in this limit while the states lc?3/2 
and lsi/2 (or I/5/2 and 1^3/2) get closer, supporting that these ones should be associated 
in the same pseudospin doublet. We notice that the assignment of the n quantum number 
to the different states is made by continuity with results at x = 0. In this case, where 
the potential is local, there is a one-to-one correspondence between this number and the 
number of nodes in the radial wave function. This correspondence is lost as soon as the 
interaction is partly non-local. Examination of radial wave functions like the lsi/2 for 
X = 1 thus evidences a node. This one occurs however in the tail of the wave function, at 
a distance where the contribution to the normalization is almost saturated (within a few 
percent's). It differs from standard nodes that allow one to determine contributions to the 
normalization from smaller and larger r of comparable sizes. Only these ones could have 
some relation to the quantum number n. We also notice that the pseudospin-symmetry 
spectrum obtained for x = 1 shows the same degeneracy pattern as the one obtained by 
Chen et al. in a relativistic approach with a harmonic-oscillator- type interaction [22] . 
Whether there is a relationship between the two cases will be discussed in Sect. 5. 



3.2 Woods-Saxon model 

For the potential V{r) inspired from the Woods-Saxon form, the two terms Vo(r) and 
Vi(r) entering its definition assume the standard form: 

Voir) = Vl(r) = -— ^ . (5) 

1 + exp((r - ro)/a)) 

For numerical calculations, we use Vq = 40 MeV, rg = 5 fm, a = 0.65 fm. Results for the 
spectrum are shown in Fig. [2l Examination of this figure confirms most of the features 
observed for the HO model. The energies of the two members of the pseudospin symmetry 
doublets, lc?3/2, 2si/2 and I/5/2, 2^3/2, are closer to each other however. As for the HO 
model, the quantum number, n, has been assigned by continuity with the spin-symmetry 
case. This quantum number is not necessarily related to the number of nodes as often 
assumed implicitly. Thus, orbits with j = I + \ tend to get an extra node when going 
gradually from the spin- to the pseudospin-symmetry limit while orbits with j = I — \ 
don't. The extra node occurs at the limit or outside the range of the potential Vi{r). This 
is illustrated in Fig. [3]for the orbits lsi/2 and 2si/2 and in Fig. IHfor the orbits 1^3/2 and 
The contribution to the norm beyond the extra node amounts to 3% in the largest 
case and tend to decrease either with going to the spin-symmetry limit or with increasing 
the orbital angular momentum. In comparison, the contribution to the norm for a usual 
node amounts to several 10% (see the two panels of Fig. [3] for a Si/2 orbit). 
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Woods-Saxon model 

Spin symmetry Pseudospin symmetry 




-0.2 0.2 0.4 0.6 0.8 



X 

Figure 2: Same as in Fig. [1], but for the Woods-Saxon model 
3.3 Comments 

Examination of the spectra for two interaction models shows that there is essentially a one- 
to-one global correspondence between states with the spin symmetry [x = 0) and those 
with the pseudospin symmetry {x = 1). There is no indication for missing pseudospin 
partners as often advocated. Instead, all states are there but in a different order. The 
order for the pseudospin-symmetry states may be a surprise but this is a consequence of 
the non locality of the corresponding interaction used here. 

The interaction model we used does not contain explicit spin-orbit interaction but it 
is found that the nonlocal term proportional to Vi provides a similar effect which, for 
X = —0.2, has the appropriate size. This is not unexpected as there is some similarity 
of this term with the one that produces the effect in relativistic mean-field approaches 
where it is proportional to a-p Vi{r) a-p. We also notice that the actual description of 
nuclei {x = —0.2) is closer to the PS limit {x = 0) than to the PSS limit (x = 1.0) and, 
moreover, requires an opposite-sign value of x. 

Our last comment is motivated by the comparison of the two models concerning the 
standard pseudospin doublets. While results for the WS model at x = —0.2 supports 
such doublets, this is not the case for the HO model. It is known from the Nilsson model 
that a term proportional to —n fi could be necessary to improve the single-particle 
spectrum, lowering, for instance, the energy of the 1^3/2, 1(^5/2 states with respect to the 
2si/2 one. This feature shows that the form of the potential is important in getting the 



7 



2s orbit 




Figure 3: Appearance of an extra node when going from the spin- to the pseudospin- 
symmetry hmit for orbits lsi/2 and 2si/2 (left and right panels respectively) with the 
Woods-Saxon model. It is noticed that the extra node appears in the tail of the wave 
function and cannot be confused with the usual nodes that appear at smaller distances, 
as the comparison of the two panels can show. The wave functions are normalized as 
Jdrp{r) = l. 



standard pseudospin doublets. This contrasts with the explanation of these doublets in 
the relativistic mean- field approach [6], which assumes that the approximate PS symmetry 
observed in this approach is, essentially, a direct consequence of the small size of + Sq 
(in comparison with the magnitude of — Sq [12]), independently of its form. 



4 Results and discussion for the wave functions 



In discussing pseudospin doublets, examination of the corresponding wave functions is 
thought to provide an important information on the relevance of the possible underlying 
symmetry. In this section, we present results based on wave functions for the standard 
pseudospin-symmetry partners as well as those resulting from our interpretation of the 
doublets. This is done in both momentum and configuration space, for the two interaction 
models we considered and for the two doublets that their spectra shown in Figs. [T] and [2] 
exhibit. 



4.1 Results in momentum space 

The transformation of the usual spin space to the pseudospin one, given by the operator 
(T ■ p, is local in the momentum space, where its effect can therefore be more easily 
calculated. Thus, the spin-orbital angular momentum part of wave functions for the 
states (i3/2, Si/2 (or P3/2) transform into that for the states P3/2, Pi/2 (or c?5/2, c?3/2), 
leaving unaffected the dependence on the modulus of p, which we denote f{p). In the 
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pseudospin-symmetry limit, these functions for the states P3/2, pi/2 (or ^5/2, ^3/2) are 
the same. Hence, in this hmit, the functions f{p) for the states ^3/2, Si/2 (or /5/2, P3/2) 
should be equal. Moreover, the functions for the states ^3/2, pi/2 (or ^5/2, ^3/2) can be 
calculated in this limit. For our interaction model, it is simply given by the functions for 
the states P3/2, P1/2 (or ^5/2, ^3/2) in the spin-symmetry limit. We thus get relations like 
the following ones: 

hd;/.XP) = hs,/.XP) = /lpi/2,3/2(P) > 

/1/5/2 (P) = /lP3/2 (P) = /w3/2,5/2 (P) • (6) 

One can therefore compare the f{p) functions for the pseudospin doublet, according to 
our interpretation, lsi/2 (or I/5/2, IP3/2), for the spin doublet 1^3/2, lpi/2 (or 

1^3/2), and for the standard pseudospin doublet l(i3/2, 2si/2 (or I/5/2, 2^3/2) (where 
we substituted the state 2si/2 (or 2^3/2) for the state lsi/2 (or 1^3/2)). For the purpose 
of the comparison, the functions for the pseudospin doublet c/3/2, Si/2 (or /5/2, P3/2) are 
calculated at x = 0. They fulfill the following normalization condition: 

Jdpf{p) = l. (7) 

The corresponding results are shown in Figs. [5] and [6] for the lowest and next to the lowest 
pseudospin doublet respectively. In each case, results for the two interaction models are 
shown. 

A quick look at the results with the two interaction models, as well as for the two pseu- 
dospin doublets, according to our definition, does not evidence much qualitative difference 
in the range which mainly contributes to the normalization. As the pseudospin symmetry 
is not fulfilled at the point x = where the functions P3/2, Pi/2 (or 6/5/2, ^3/2) are calcu- 
lated, one expects that the various curves shown in the figures evidence some discrepancy. 
With this respect, these discrepancies are not small. At low p, they essentially reflect the 
p^~^^ behavior expected for a local interaction model. At higher p, they rather reflect the 
average momentum of the state, which increases with its energy as measured from the 
bottom of the potential well. It is nevertheless seen that functions f{p) for the pseudospin 
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HO model 



WS model 




P [fm '] p [fm '] 



Figure 5: Comparison of momentum-space wave functions, f{p), for the lowest pseudospin 
doublet with the HO and WS models (left and right panels respectively): small- and 
large-dash lines for the states lsi/2, 1(^3/2, continuous line for the limit expected in the 
pseudospin-symmetry limit and dotted line for the 2si/2 state. All wave functions shown 
in the figure are calculated at x = 0. 

doublets, as considered here, have a unique sign in the range which mainly contributes to 
the norm. Moreover, they share this feature with the function expected in the pseudospin 
symmetry limit. Instead, the functions for the standard pseudospin doublets evidence a 
different structure in the same range. One of the function has a given sign while the other 
one (with a different n number) evidences a change in sign. It thus turns out that our 
assignment for pseudospin doublets is more in agreement with the pseudospin-symmetry 
expectation than the standard assignment. 



4.2 Results in configuration space 

While the direct comparison of wave functions in momentum space makes sense, it does 
not in configuration space where wave functions imply an integral over p of momentum 
wave functions together with spherical Bessel functions of different /. A possible compar- 
ison could involve the integral of the different momentum wave functions with the same 
spherical Bessel function ji(pr) for the lowest pseudospin doublet (or j2{pr) for the second 
one). We defined these quantities as: 

/(^) = \/|- ^ / f(P) 3i{pr) (or j2{pr)) , (8) 

which, using Eq. ([7j), are found to verify the normalization condition: 

jdrf\r) = l. (9) 

Interestingly, the quantities /(r) together with the associated angular momentum struc- 
ture corresponding to the states P3/2, Pi/2 (or ^5/2, (^3/2), show some similarity with the 
small components in the Dirac mean-field phenomenology. It is reminded that the ex- 
amination of the equation allowing one to determine these components provided some 
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HO model 



WS model 




Figure 6: Comparison of momentum-space wave functions, f{p), for the second lowest 
pseudospin doublet with the HO and WS models (left and right panels respectively): 
small- and large-dash lines for the states 1^3/2, 1/5/2, continuous line for the curve ex- 
pected in the pseudospin-symmetry limit and dotted line for the 2^3/2 state. All wave 
functions shown in the figure are calculated at x = 0. 



justification for the existence of pseudospin doublets [6l[TT]. Taking into account Eq. 
the functions /(r) are expected to verify the following relation in the pseudospin-symmetry 
limit: 

/1/5/2 (r) = /ip3/, (r) = /id3/2,5/2 (r) ■ (10) 

Results showing the /(r) functions are presented in Figs. [7] and [HI in full correspondence 
with those shown in Figs. [5]and[6]for momentum space. Qualitatively, the different curves 
evidence a similar pattern. There is not much sensitivity to the interaction model and the 
slight sensitivity to the pseudospin doublet can be mainly ascribed to the corresponding 
binding, the results for the less bound system tending to extend to larger radial distances. 

In each panel, curves corresponding to the pseudospin assignment made here have the 
same sign within the largest range of their contribution to the normalization given by Eq. 
(Q . They evidence some spreading but fall slightly apart on one side and on the other side 
of the expected result in the pseudospin-symmetry limit. The curves corresponding to the 
standard pseudospin assignment evidence a striking difference. One of them changes sign 
within the range of its main contribution to the normalization. Thus, these results show 
that, like in momentum space, if one considers only the similarity of the /(r) functions of 
the PS partners, it is more natural to consider the states 1^3/2, lsi/2 (or I/5/2, 1^3/2), as 
a pseudospin doublet rather than the states 1^3/2, 2si/2 (or I/5/2, 2^3/2). The argument 
will get further strength if we notice that the above spreading represents for a large part a 
difference in the binding energy or in the momentum distribution of the wave functions. It 
is found that Eq. (fTOj) is much better fulfilled than what Figs. [7] and [8] suggest if we require 
that the average square momentum be the same, as expected in the pseudospin-symmetry 
limit. This can be checked by using scaled wave functions. 
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HO model WS model 




6 8 10 2 4 

r [fm] r [fm] 



Figure 7: Comparison of configuration-space wave functions for the lowest pseudospin 
doublet with the HO and WS models (left and right panels respectively). See caption of 
Fig. Ofor the definition of the different curves. 



HO model WS model 




r [fm] r [fm] 

Figure 8: Comparison of configuration-space wave functions for the second lowest pseu- 
dospin doublet with the HO and WS models (left and right panels respectively). See 
caption of Fig. [H] for the definition of the different curves. 

5 Relation with results from a relativistic approach 

Since the PS symmetry has been considered in the last decade as a symmetry with deep 
roots in the relativistic theory of atomic nuclei [6j, arises the question of whether the 
above results can be cast into this approach, which can be described by the following 
equations: 

a-p ^<(f) -{E - M - V{r)) V^>(f) = 0, 

(T-p ^>(f) - (^ + M- A(r)) ^<(f) = 0. (11) 

where ipy{r) and ip<:{r) represent, respectively, the usual large and small component of 
the Dirac spinor. Various models are considered in the following. 
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5.1 Model inspired from the non-relativistic one 

Generally, the potentials V{r) and A(r) are taken as the sum and the difference of the 
scalar and vector potentials, V{r) = V^{r) + Vs{r) and A(r) = V^^r) — Vs{r) , but , formally, 
nothing prevents one to take V{r) as given by Eq. ([3]) and A(r) as a constant. This model 
has not much to do with the usual relativistic mean-field phenomenology but it can be 
easily studied. Actually, it does not add relevant information to the non-relativistic model 
studied in the previous section. Very similar results are obtained for the spectrum as well 
as for the wave function as far as the large component is concerned. 



5.2 Model inspired from the Chen et al. one 

A more interesting comparison involves the model considered by Chen et al. [221 123] . 
Apart from the energy scale, the pseudospin-symmetry spectrum shown in Fig. [1] evi- 
dences a striking similarity with theirs. There could be also some similarity with results 
obtained by Lisboa et al. [36] , which are essentially the same as those obtained by Chen et 
al., but the authors modified the assignment of the quantum number n to make it identi- 
cal to the number of nodes of the radial wave functions. The authors compare the results 
in the cases of the two symmetry limits, arguing about missing and intruder states. The 
question therefore arises whether there is a continuity between the two limits, as observed 
in the non-relativistic case discussed in previous sections. Such a study could also help in 
clarifying the assignment of the n quantum number. 

In order to study the relationship between the two symmetry cases, one has to determine 
the two potentials, V{r) and A(r), that enter the following standard equations for the 
large and small components, G{r) and F{r): 



^K'±ll-(^E+M-A{r)){E~M-V{r)) 



lihi) 



- {E+M-A{r)) {E-M-V{r)) 



G{r) 
F{r) 



E + M-A{r) 

dr ^ ar r ' \ / / 

E-M-V{r) 



0, (12) 
0, (13) 



where the notation such as (;f; + 7, G{r)) stands for (G'(r) + ^ G{r)). In the cases of spin- 



and pseudospin-symmetry limits, the potentials V{r) 
given by: 



and A(r) should be respectively 



V{r) 



Voir) 
V{r) 



A(r) = 0, 



1 



A(r) = Ao(r) = ^tu^Mr^ 



(14) 



A first guess to determine the transition potential would consist in multiplying Voir) 
and Ao(r) by the factors \ — x and x, and allow x to vary between and 1. It can 
be however shown that solutions of Eqs. (fT2| [T3!) have then an oscillatory behavior at 
large r, which complicates their determination and could obscure conclusions. As we 
are interested in finding at least one way to ensure a continuous transition between the 
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two symmetry limits, we proceed differently. Instead of fixing a priori the transition 
potential, we determine this potential from minimal consistency conditions. To illustrate 
the method we adopted, we consider here with some detail the case n = 1 and 1 = 1 + 1. 

First, as the above mentioned oscillatory character of solutions is due to the appearance 
of an attractive term in the product {M+E — A{r)) (M — i?+V^(r)), we skip this term 
so that the remaining part behaves as an ordinary harmonic-oscillator potential. This 
product thus read: 

{E+M-^{r)){E-M-V{r)) = E^-M^ - ^ ((l-x)(E+M) + x{E-M)) . (15) 

It allows one to determine V{r) once A(r) has been determined from other sources. 
Second, the examination of the last term in Eq. (|T2|) in the pseudospin-symmetry limit 
{V{r) = 0) shows that the singular behavior due to its denominator is exactly cancelled 
by the action on G{r) of the operator, ^ + 7, at the numerator. The equation can then 
be solved analytically, allowing one to recover the solution which is known in any case. 
It is not difficult to see that the procedure can be extended to any x in the simplest case 
of orbits with n = 1 and / = / + 1 (states lpi/2, 1(^3/2, I/5/2, l5'7/2, ■ • ■)• The solution 
assumes the form: 

G(r) oc r('+^) exp(-«VV2) , (16) 
while the last term in Eq. (fT2l) can be written as: 



E + M-A{r) {E+M-A{r))r E+M 

By inserting the above expression in Eq. (fT2|) . one gets the following relations to be 
fulfilled: 



a"^ -cui —({l-x){E + M) + x{E-M)) =0, 



2 



E^-M^ - u (2/ + 3) ((l-x)(E+M) + x{E-M)) + + = . (19) 

The first of them is obtained from cancelling the coefficient of the highest r-power term in 
the equation determining G{r) (term proportional to r'^^). The second one allows one to 
get the energy and is obtained from cancelling the next term in the same equation (term 
proportional to r'"*"^). It can be checked that the states lpi/2, 1(^3/2, I/5/2, ^97/2, • ■ • in 
the spin-symmetry limit are in a one-to-one correspondence with the same states in the 
pseudospin-symmetry limit. Assigning to these states the same n quantum number as 
for the spin-symmetry case, it is found that this number coincides at x = 1 with the one 
assigned by Chen et ai, in this case from algebraic considerations [22] . 

For completeness, we give the potentials A(r) and V{r). The potential A(r) is obtained 
by solving the equation that results from the last equality in Eq. ([1] 



^ Lj^M (21 + 1) r , , 

E+M-A{r) E+M {2l + l-a^r^)' ^ ' 
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The solution of this equation is given by: 



„.2 2 x(2l+l)<^'^M 

E+M-A{r) = (E+M) (l-^^) • (21) 

In integrating Eq. fl20|) . we fixed the integration constant by requiring that A(r) = 0, 
as expected in the spin-symmetry hmit corresponding to a; = 0, hence the overall factor, 
E + M, at the r.h.s. of Eq. fl2T|) . The other factor is well defined when the quantity, 
2/ + 1 — a^r^, is positive. For negative values, it is appropriate to replace this factor by 
e^*'^(Q;^r^ — 2/ — 1). The imaginary phase that could then occur, exp(±27r^|^^^^), can 
be removed but should be accounted for in determining the small component. It can be 
checked that the exponent is equal to 1 at x = 1, so that the factor at the denominator 
in Eq. f|T7j) is exactly cancelled by the factor (2Z + 1 — a^r^) at the numerator (up to a 
constant factor). The expression of the small component is easily obtained and reads: 

= (f + 7> Gj-)) ^ m (22) 

^ ' E + M-A{r) E+M V 2/+l^ r ' ^ ' 

The potential V{r) can be obtained from the expression of A(r), given by Eq. fl21l) . and 
from Eq. (fT5l) . which can be cast into the form: 

{M+E-A{r)) {M-E+V{r)) = M'^-E^ + a^r^. (23) 

The above approach can be extended to other orbits, either with non-zero n or with 
/ = /-|-1. We give below equations allowing one to determine as a function of x and 
make some comments about obtaining these results. We refer to the appendix |A] for 
further details. 

In all cases, we assume that the wave function writes as the product of an exponential 
factor, exp(— Q;^r^/2), and a polynomial term, similarly to the harmonic oscillator wave 
functions that are recovered at x = for the large component, G{r), or at a; = 1 for the 
small one, F{r). The consideration of the highest r-power term in Eq. 0121) ensures that 
the expression of has always the form given by Eq. f|T8|) . while the consideration of the 
lowest r-power term in the equation requires that the wave function contains the minimal 
well known factor r'+^. 

For 1 = 1 — 1 {k = I) and any n {npi/2, ndz/2, ngj/2, ■ ■ ■), we assume that the 

polynomial part of the wave function has a structure similar to the spin-symmetry case 
and, in particular, has the same number of nodes. From considering the next to the 
highest r-power term (r^*^"^) in Eq. ( fT2i) . one gets the equation giving E' as a function 
of x: 

E^-M^ - {4n+2l-l) + x- — ^ = 0. (24) 

E + M 

The lowest r-power terms are useful in determining the polynomial part and, from it, the 
expression of the potential A(r). This is tractable for the case n = 1 considered above 
(first-degree equation) and for the case n = 2 detailed in the appendix |A] (second- degree 
equation) but the difficulty of the task quickly increases with n (n-degree equation). 
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0.2 0.4 0.6 0.8 



X 

Figure 9: Single-particle spectrum in a relativistic approach with a harmonic-oscillator 
inspired interaction: transition from the spin-symmetry limit (x = 0) to the pseudospin- 
symmetry limit (x = 1). We limited the representation of states to those corresponding 
to the complete multiplets shown on the r.h.s. of the figure. Some higher-energy states 
are therefore not shown. 



For finding solutions for the case with / = / + !(/? = — (/ + 1)) and any n 
{nsi/2, np3/2, nd^i2, nfT/2, ■ ■ ■), we notice that there is some symmetry between Eq. 
f[T^ and Eq. f[T^ . The orbital angular momentum / is replaced by / and the solution for 
the small component has the same number of nodes at x = and x = 1. What we did 
previously for the large component can therefore be repeated here for the small one with 
the appropriate replacements (A(r) V{r), M —M, ). Thus Eqs. 

(120]) and ([21]) for n = 1 become: 



dV{r) 
dr 



(!^, (25) 



E-M-V{r) ' ' E-M {2l + l-a^r^) 

and: 



E-M-V{r) = {E-M) [l-^^--^) s.^c^^m) _ (^26) 
The general equation allowing one to determine the spectrum for any n now reads: 

E^~M^ -{^n\2l-\) \{\-x)- — —- = 0. (27) 

E — NL 
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The results for the spectrum are shown in Fig. [Hlfor x varying from to 1. They have 
been obtained using parameters identical to those used in Refs. ([23], [36]) (M = 10 fm~^ 
and uj = 2 fm~^), so that to facilitate some comparison with their results. This choice 
also offers the advantage that the spectrum at a; = 1 is less compressed as for the choice 
uo = 1 ImT^ used in an earlier work [22] . Examination of the figure shows that there is 
an one-to-one correspondence between states in the spin- and the pseudospin-symmetry 
limits, indicating that there is no missing states as one goes from some limit to the other 
one. 

lS|^, orbit Ip^^^ orbit 




X [fm] r [fm] 

Figure 10: Appearance of an extra node in the component G{r) when going from the 
spin- to the pseudospin-symmetry limit for orbits lsi/2 and lp3/2 (left and right pan- 
els respectively) with the Chen et al. model. The normalization assumes the relation 
!dr {G^{r) + F^{r)) = 1. 



In relation with the wave functions, we recover their expressions expected in the two 
symmetry limits. We find, in particular, that the large component of the states with 
/ = / + 1, corresponding to a given n (as defined here by continuity with the one for 
the spin-symmetry case), acquires an extra node while one goes from the spin- to the 
pseudospin-symmetry limit. It shows that establishing a correspondence of the states with 
the same number of nodes in both limits, as done in Ref. [36], could lead to misleading 
conclusions. Actually, examination of the wave functions shows that the extra node 
appears at a distance where the contribution to the norm is almost saturated. This result 
is similar to the one obtained in the previous section for a non-relativistic model but the 
rate of saturation is smaller in the present case. Contributions to the norm beyond the 
last node amount to around 10% for the lsi/2 state, 7% for the lp3/2 state and less for 
higher angular momenta. Some results are shown in Fig. (TDlfor the lsi/2 and lp3/2 states 
while expressions for the wave functions can be found in the appendix [XI Contrary to 
the nonrelativistic case, the appearance of an extra node is not progressive. Virtually, the 
node is present at a; = but it does not show up because the factor responsible for it, 

2 2 

1 — is taken with the power 0, making this factor equal to 1 for any r value. This 
is a particular feature of the transition model we used. The appearance of an extra node 
also holds for the small component, but only for the states with / = l—l and when making 
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the inverse path. 



We assumed in the above study of the continuity between the spin- and pseudospin 
symmetry hmits that the last term in Eq. ( |T2l) was hnear in x, partly because it was the 
simplest possible assumption to start with. Many other assumptions could be made. A 
particularly attractive one for the case I = l—l supposes to replace the factor —x '-^''^^^^ 

at the r.h.s. of Eq. (ITTI) by — 2x^^^^^^. The dependence on x of this expression is 
more complicated than for the choice retained here {a also depends on x) , but the result 
for the relation between E and x turns out to be formally simpler. It is given by: 

E^-M^ = {'in+2l-l-2x)a\ (28) 

In writing the above result, we anticipated its generalization to any n. A similar result 
could be obtained for the case / = /+1 with the replacements of Z by Z and x hj 1 — x. It 
reads: 

E^-M^ = {4:n+2l-l + 2x)a'^ . (29) 

The continuous transition between the spin-and pseudospin-symmetry limits is somewhat 
straightforward from the above expressions. In both cases, the result should be completed 
by the derivation of the potentials V{r) and A(r), which is not necessary simpler than 
before. 



5.3 Model with V{r) and A(r) inspired from Woods- Saxon ones 

In previous sections, we examine models that were providing some continuity between 
the spin- and pseudospin-symmetry limits. There is no guarantee that such a continuity 
exists in all cases however. On the one hand, some states may disappear and reappear 
(or appear and disappear) in going from one limit to the other but in cases considered up 
to now we could find a model avoiding this drawback. On the other hand, some states 
may disappear (or appear), what cannot be excluded in a relativistic model where the 
spectrum is not bounded from below. A reason to be suspicious comes from the fact that 
realizing the pseudospin symmetry with a Woods-Saxon type potential supposes that its 
depth is larger than twice the nucleon mass, which is the obvious sign of a relativistic 
regime. 

We examined various transition models from the spin-to the pseudospin-symmetry lim- 
its but, due to various problems (numerics, interpretation), we could not get clear con- 
clusions in most cases. We present here results for one case where some interpretation 
is possible. Moreover, they provide interesting clues as for what could occur between 
the two symmetry limits. The model is described in the appendix while the spectrum is 
shown in Fig. [11] An essential difference with Chen et aZ.'s model is that the potentials 
at X = and x = 1, instead of going to oo with r, are now limited by an upper finite 
value. 

We first notice that the spectrum in the spin-symmetry limit looks very much like the 
usual one. As far as the order of the states is concerned, the spectrum in the pseudospin- 
symmetry limit is similar to the one for the WS model shown in Fig. [2J It would also 
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Figure 11: Single-particle spectrum in a relativistic approach with a Woods-Saxon inspired 
interaction: transition from the spin-symmetry limit (x = 0) to the pseudospin-symmetry 
limit (x = 1). Notice that our conventions for the potentials at x = and x = 1 imply 
positive values for the quantity, E — M, represented on the vertical axis (see appendix 

be similar to the one obtained by Chen et al. [23], had the authors chosen to raise the 
degeneracy of their HO interaction model with a Woods-Saxon type correction having 
the opposite sign. Considering the large component of the wave functions, it is found 
that all states with j = 1 — 1/2 and same n quantum number (defined here by the 
order of appearance in the spectrum) have the same number of nodes in the pseudospin- 
symmetry limit as in the spin-symmetry limit. Instead, states with j = / + 1/2 get an 
extra node, located however in a domain where the contribution to the norm beyond this 
node is, generally, small. This observation is similar to the one we could make for the 
non- relativistic case or for Chen et a/.'s model [22] . 

At first sight, from the qualitative agreement of the spectra in the SS and PSS limits 
with those obtained previously in this work, we could expect similar conclusions for the 
spectrum between these two limits. Looking at this one however reveals important dif- 
ferences. It is found that there is continuity between x = and x = 1 for all states with 
j = 1 + 1/2, without change for the n quantum number (always defined here by the order 
of appearance in the spectrum). This is also satisfied for states with j = / — 1/2 and / > 3 
but the picture differs for states npi/2 and As these states evolve from the spin- 

to the pseudospin-symmetry limit, their n quantum number increases by one unit. This 
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change is related to the appearance of an extra node in the core of the large component 
of the wave function. As a counterpart of the change in n, the states lpi/2 and 1^3/2 in 
the pseudospin-symmetry limit cannot be related to any state in the spin-symmetry limit. 
As seen in the figure, they disappear from the spectrum around x — 0.5 — 0.6, where the 
interaction is strongly modified. This fact explains the disappearance of the PS partner 
of the intruder lsi/2 state. It is also seen on the figure that a new kind of states with 
j = 1 + 1/2 (to whom we assign, arbitrarily, n = 0) appears and disappears between ,x = 
and x = 1, without spoiling the symmetry properties of the spectra obtained in these 
limits. 

One can think, perhaps, that the extra states that appear with > for values of 
a; > are simple artefacts of the interaction model. However, among them, those states 
with j = 1 — 1/2 are necessary in order the PS symmetry be satisfied for x = 1. The 
other states with j = I + 1/2 are not required by this symmetry but their appearance 
strongly influences the spectra of the states that arc above, pushing them upwards and 
evidencing level-repulsion properties around x = 0.5 — 0.7. Level-repulsion properties are 
also observed for the levels with j = Z — 1/2 in the same range. In these cases, this is 
the range where the wave function can be strongly changed and acquires an extra node. 
In some of these cases, the small and large components become comparable, suggesting 
the role of relativity in this change. It is noticed that, if we ignore the small region 
where the level-repulsion eff'ects for these states appear, the spectra look continuous from 
the spin- to the pseudospin-symmetry limits. Properties of the large component of the 
wave function before and after the pseudo-crossing point with a similar energy are not 
much affected. These features are similar to the ones for the dependence of single-particle 
spectra as a function of the deformation of the nucleus. 

The detail of the above results depends on the transition model but part of it, hke the 
appearance of an extra node as x increases as well as the appearance (or disappearance) 
of some states from the positive part of the SP spectrum, may point to relevant features 
of the models used to describe the interaction in the spin- or pseudospin-symmetry limit. 
Further studies could be useful to determine, in particular, whether the disappearance of 
states could be completely removed or, on the contrary, be extended to all states with 
j = I — 1/2 and n = 1, then providing support for the missing partners of intruder states. 
Whatever the picture, the results point to an important breakdown of the pseudospin 
symmetry. States disappear from the spectrum in one case or, in the other case, the 
assignment of states to pseudospin doublets differs from the one usually assumed, while 
agreeing with results from other parts of this work. 



6 Conclusions 

In this paper, we examined whether there could be some continuous transition between 
the spin-symmetry and the pseudospin-symmetry limit description of atomic nuclei and 
looked at what would be the implications of such a result. This was done by introducing 
an interaction that is a continuously varying superposition of potentials satisfying these 
symmetry properties. 
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Within a non-relativistic approach, we found that there could be a one-to-one corre- 
spondence between the spectra calculated in the SS and PSS limits. This is obtained 
with an interaction that is non-local in the PSS limit, resulting in a state ordering of the 
spectrum different from the one usually expected (the ground state is not necessarily a 
Si/2 one). For the transition model we considered, energy shifts are larger than for the 
spin-orbit splitting but they do not generally exceed the splitting corresponding to two 
major shells. The results thus suggest to associate in the same PS doublet states lsi/2 
and IP3/2 and I/5/2, ■ ■ ■, rather than states 2si/2 and 1^3/2, 2^3/2 and I/5/2, • ■ ■• Our 

suggestion is further supported by the analysis of wave functions which leads to results 
closer to the PS-limit expectations with our assignment than with the usual one. With 
this respect, we notice that wave functions for the states lsi/2, IP3/2, ■ • ■, get an extra 
node in the PSS limit but this node is located in the tail of the wave function rather 
than in its bulk as for states 2si/2, 2^3/2, • ■ •. This prevents one to identify these states 
together, what is further supported by our study looking at a continuous transition from 
the SS to the PSS limit. 

We also looked at a possible continuous transition between the spectra calculated in the 
SS and PSS limits within a relativistic approach as this one was thought to provide some 
support for PSS doublets. Without regard to the origin of the interaction, we first notice 
that previous results can be cast into the relativistic approach by identifying the potential 
V{r) to the one used in the non-relativistic case and taking A(r) as a constant. Another 
relativistic model is the one considered by Chen et al. p2j, based on oscillator type 
potentials. In this case, we were able to find a potential allowing one to get a continuous 
transition between the SS and PSS spectra, which turn out to be very similar to those we 
get in the non-relativistic case. Contrary to statements made by the authors, there is no 
need to include intruder or missing states to make the spectra consistent with each other. 
We simply have to accept that the order of the states in the PSS limit differs from the 
usual one. With the intent to study more realistic models, we finally consider the case of 
finite-size potentials. Apart from the energy scale, the spectra in the SS and PSS limits 
were found to be similar to those obtained in the previous cases, suggesting that one could 
establish some continuity between them as done previously. However, while all states in 
the SS limit could be continuously related to states in the PSS limit, the reciprocal was 
not true. Some states in the PSS limit could not be related to states in the SS limit and, 
actually, were found to disappear from the positive-energy spectrum. Interestingly, the 
spectrum so obtained is found to be similar for a part to the observed one, including the 
assignment for PSS doublets. Thus, these results could throw light, in particular, about 
the disappearance of the PS partners of the intruder states when one goes from the PSS 
to the SS limit. The results nevertheless call for two important remarks. First, as the 
disappearance implies a limited number of states with n = 1 and j = / — 1/2, there is 
likely a dependence of the results on the interaction models we used in the SS and PSS 
limits. Second, the disappearance of states from the positive-energy spectrum points to a 
very strong violation of the PS symmetry in this case, more important than the one found 
in the other approaches where only a re-ordering of the levels was observed. Moreover, the 
assignment for PS doublets then implies that the large component of the wave function 
of one of the partners gets an extra node in its bulk when going from the SS to the PSS 
limit. This feature represents another violation of the PS symmetry. 
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From examining wave functions, it was shown in the hterature that the pseudospin 
symmetry could not be as good as the spin one [281 EOl [32], pointing to the accidental 
character of the validity of this symmetry to describe the SP-level spectrum. Starting from 
a different viewpoint, the present results largely support this observation. If some doubt 
is allowed, it is due to one of the studies we performed within a relativistic framework. In 
this case, however, further studies could be useful to determine in what extent the results 
depend on specific features of the interaction we used. Either these results are confirmed 
and they imply, intrinsically, a very strong violation of the pseudospin symmetry, or, they 
confirm results of the other approaches, with a violation of this symmetry that is smaller 
but nevertheless larger than the spin- symmetry one. 
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A Transition from the SS to PSS limit in the Chen 
et al. model 

We give in this appendix some further details relative to the calculation of the spectrum 
between the spin- and pseudospin-symmetry limits in the Chen et al. model [22] • They 
successively concern the equations relating E to x, and the parameters u and M, the 
determination of the transition potential for the case n = 2 and / = l — l, and, somewhat 
sketched, for the case n = 1 and 1 = 1 + 1. 

• Remark about the relation of E to u 

It is first noticed that the equations allowing one to determine the relation of E to x, 
Eqs. (l24l) and (|27l) . are quadratic in u and, therefore, apparently differ from those given 
in Ref. p2] in the case x = 1 and x = respectively, which are linear in u. They can 
nevertheless be rearranged so that to express u in terms of E and x. When this is done, 
one recovers the expressions given in the above work: 



u(2n + l-^)=u{2n + l-^) = {E+M)J^^ {1 = 1-1, x = l), (30) 



u;{2n + l-^)=u;{2n + l-^) = iE-M)J^±^ {1 = 1 + 1, x = 0). (31) 



• Case of states with n = 2 and j = I — 1/2: 2pi/2, 2^3/2, 2/5/2, ■ • (1 = 1 — 1) 

We assumed that the wave function keeps the same form as for a; = or x = 1, G{r) oc 

r'^'+^^(a^r^ — y) exp(— r^/2), which implies that the last term in Eq. ([T2D can be written 

as: 



^(f + f,G(r)) ^ (a^T\2l + ?>+y)-{2l + l)y-a'T' 



E + M-A{r) {E + M-A{r)){a'^r^ -y)r 

(2/ + 5)w2M (aV^ 



G(r) 



-x^— ^— — f— ^ zGir). (32) 

E+M {a^r^ -y) ^ ' ^ ' 
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The last equality is a further assumption which implies an equation to be solved: 

^ ^ (21 + 5) u'M rja^r'-z) ^^^^ 



E+M-A{r) E+M (^2 ^2(2/ + 3+y) - (2/+ 1) y 



a4 r4 



The quantities y and z are two constants that have to be determined. The first one, which 
corresponds to a node in the large component of the wave function, is known at x = 
{y = and at X = 1 {y = ^^). The second one corresponds to a node in the small 
component of the wave function and is known at x = 1 where it has a simple expression 
{z = ^^)- It could be chosen to simplify part of the calculations. Both y and z enter 
Eq. (fT2|) to be solved. Looking at this equation, it is found that the coefficients of the 
terms proportional to r'"*"^ and r'~^ vanish with the choice we made for the wave function 
(using Eq. (fT8|) in the first case and the low-r behavior, r'"*"^, in the second one). The 
vanishing of coefficients of terms r'+^ and r'"*"^ supposes the conditions: 

9 9 . , ^ 9 (21 + 5) U^M 

E^-M' - (21 + 7) a' + x ^ = , 

(E^-M^) y - (2/ + 3)(2 + y) + x z = 0. (34) 

The first of them can be cast in the general form given by Eq. ([23]). The second one allows 
one to get y once a choice has been made for z. This last quantity enters the determination 
of A(r) and, therefore, the potential which allows one to make the transition between 
X = and x = 1, where it identifies to the one used in Chen et a/.'s work. A particular 
choice consists in taking for z the value of a'^r^ that cancels the denominator at the r.h.s. 
of Eq. (15^ . the other value being denoted z'. The value of y, z, z' so obtained are given 
by: 

y = ((2/+3) {{AEr + GAEW + 16a^) - {2l + l){^^f 



-2AEY(AE2 + (2/ + 7)a2)2 + 4(2/ + l)(AE2 + 2a2)a2 

1 



(2/ + 3)(Z^^ + 6«2) - (2/ + l)Z^2 



^a2 



z' = 2l + 3 + y-z 
1 / 



2^Z^2+4^2^ 



2 J(AE2 + (2/ + 7)a2)2 + 4(2/ + l)(AE2 + 2a2)a2 



{21 + 3){AE^ + 6a^) + {2l + l)AE' 



+2A/(AE2 + (2/ + 7)a2)2 + 4(2/ + l)(AE2 + 2a2)a2 , (35) 



where AE^ = E^-M^ - {2l + 7)a^ = -x ^ ^e+m = at x = and AE2 = -2 at 

X = 1). The quantity A(r) can then be calculated easily from the following equation: 

^ -.^^(2^ ,36) 



E + M-A{r) E+M (z'-a^r^) 
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The solution takes a relatively simple form. It is similar to Eq. fl2ip for the n = 1 case 
and is given by: 

^2 2 £(2i+5)iAM 

E+M-A{r) = (E + M) (l-^) '"''"^""^ ■ (37) 

It can be checked that the exponent is equal to 1 at x = 1, so that the factor at the 
denominator in Eq. ( l32ll exactly cancels a similar factor (2;' — a^r^) at the numerator (up 
to a constant factor). 

The solution for the small component is given by: 

pr.l = ii + V^Gir)) ^ (2m) / {l-^)Gir) 

^^^^ E+M-A{r) E+M y z'^ ' ^ > 

• Case of states with n = l and j = I + 1/2: lsi/2, 1^3/2, 1(^5/2, I/7/2, • ■ •; (1 = 1 + 1) 
The solution we retained in this case has been obtained from the observation that there 
is some symmetry between Eq. ( |T2|) for the large component and Eq. ( |T3l) for the small 
one. Most developments presented in the main text for the case (/ = / + 1) can thus be 
transposed here. For the case with n = 1, we therefore look for a small component with 
the following form F(r) cx r'^'^-^-' exp(— r^/2) oc r^'^^-' exp(— r^/2), which is common 
to the expected solutions at x = and x = 1. The solution for the large component can 
be obtained using its expression in term of the small one together with the expression of 
the transition potential given by Eq. fl26|) : 



^ ' E-M-V{r) E-M ^ 21+V r ^ ' 

The above solution is consistent with the idea that the spin symmetry results from the 
cancellation of A(r) while the pseudospin symmetry results from the cancellation of V{r), 
without any regard to the smallness of one component with respect to the other one. 
Solutions that are more in agreement with this last observation have also been found. 
They nevertheless differ from the solutions presented in the main text with many respects. 
The equation allowing one to determine ii^ as a function of x is more complicated but the 
appearance of a node for the large component is more progressive. On the other hand, 
some modification of the normalization, that can be actually justified, is necessary to 
make it convergent. 



B Transition from the SS to PSS limit in a relativistic 
model with finite potentials 

In order to study the transition between two models A and B satisfying exact SS and 
PSS, respectively, we consider, in accordance with the definitions given in Subsect. 15.11 
the potentials V{r) = S^(r) + Ss(r) and A(r) = St,(r) — Ss(r) entering the Dirac equation 
(HH): 
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V{r) 



8M 

1 + e~(^~''o)/" 



1 



+ 1] 



1 



1 + e{r-ro)/a ' ]_ _|_ g-(r-ro)/a 



(40) 



Afr) = X- 



AM 



1 + e^('"^'"o)/'^ ' 

where, we take, somewhat arbitrarily, a = 0.1, M = 10 fm~^, tq = 2 fm, and 



4:{xo - x) 



with (see below) 



l + 100(a;-a;o)2^ 
E + M 



xo 



AM 



(41) 



(42) 



(43) 



Thus, for X = we recover the model A, which satisfies the exact spin symmetry, whereas 
for X = 1, we have the model B, which satisfies the exact pseudospin symmetry. We do not 
think our main conclusions will depend very much on our specific election of parameters. 
Like for Chen et al.'s model [22], the strengths of the potential V{r) at x = and A(r) 
at X = 1 are small in the vicinity of r = 0. In the large-r limit, they are however bounded 
by an upper finite value instead to increase infinitely. The value of e = E' — M is then 
necessarily positive at x = or x = 1. In what follows, we may omit, by simplicity, the 
explicit dependence of V{r) and A(r) on the space coordinate r. 

To facilitate the understanding of the solutions of the Dirac equation (ITT]) with these 
potentials, it is useful to consider the corresponding equivalent Schrodinger equation for 
the large component (see Ref. [32] for details). In this equation, with the potential A 
that becomes constant for r — > oo, the central potential {U{r) = t/) at r — > oo can be 
approximated as 



f/(r 



oo 



(44) 



where e = E — M. 



The last term in Eq. fHOj) . oc rj, makes V{r — > oo) different from zero. Its sign changes 
at X = Xq to make the central potential U{r) attractive at large values of r. The factor 
4 entering the expression of r] controls the "speed" of the transition from the SS to the 
PSS limit as x varies. The factor 100 (together with the factor 4) determines the value 
of at r — i> oo for a given value of x. 



The potential V{r oo) can be written in terms of U and A as, 

V{r 



, 2M{U/e) - (A - e) 
oo) = — — — : e. 



2M - (A - e) 



(45) 



Notice that Xq is the value of x for which the denominator 2M — (A — e) = 0. The 
potential V is chosen so that it remains always finite at r — > oo, in particular for x = Xq. 
In this case, necessarily, U{r oo) = e for Thus, the system for X — ^ Xq spreads 
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all over the space. However, with the choice for the potentials V and A, it is possible 
to find bound states for x very close to xq and, in most of cases, the SP energies vary 
smoothly with x around xq- Only when the value of the SP-state energy is high enough 
and Xq approaches the true value of x, for which the transition from the SS to the PSS 
limit is produced, in practice, the SP energy may vary quite sharply for states with / large 
and j = I — 1/2 (see the SP energies for the states /5/2 and gj/2 in Fig. [TTi) . 
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